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ABSTRACT 

We consider the evolution of white dwarfs with compact object companions (specifi- 
cally black holes with masses up to ~ 10^ Mq, neutron stars, and other white dwarfs). 
We suppose that the orbits are initially quite elliptical and then shrink and circularise 
under the action of gravitational radiation. During this evolution, the white dwarfs 
will pass through resonances when harmonics of the orbital frequency match the stellar 
oscillation eigenfrequencies. As a star passes through these resonances, the associated 
modes will be excited and can be driven to amplitudes that are so large that there is a 
back reaction on the orbit which, in turn, limits the growth of the modes. A formalism 
is presented for describing this dynamical interaction for a non-rotating star in the lin- 
ear approximation when the orbit can be treated as non-relativistic. A semi-analytical 
expression is found for computing the resonant energy transfer as a function of stellar 
and orbital parameters for the regime where back reaction may be neglected. This is 
used to calculate the results of passage through a sequence of resonances for several 
hypothetical systems. It is found that the amplitude of the £ = m = 2 /-mode can 
be driven into the non-linear regime for appropriate initial conditions. We also dis- 
cuss where the no back reaction approximation is expected to fail, and the qualitative 
effects of back reaction. 
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1 INTRODUCTION 

White dwarfs, the normal evolutionary endpoint for stars 
less massive than ~ 8 M0, are extremely common; the halo 
of our Galaxy contains several billion of them. They are 
observed frequently in binary systems, with normal stellar 
companions, as cataclysmic variable stars and, less often, 
with compact object companions, WDCO systems. Many of 
these systems are produced naturally in binary star evolu- 
tion and, as a consequence, have circular orbits. However, it 
is also possible to form eccentric, WDCO binaries following 
stellar capture or exchange in a dense stellar environment. 
In the present paper, we do not consider the formation of 
such binaries in detail; we list some possible formation mech- 
anisms, but our detailed considerations are restricted to the 
evolution of such systems after they form. 

There are three possible types of WDCO systems. 
Perhaps the most common are double degenerates (e.g. 
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iMaxted et alJl2002f) . which have been observed to have pe- 
riods as short as ~ 5 min ([Ramsav ct al. 2002). A second 
class is white dwarf-neutron star binaries, of w hich there 
are over 50 known or suspec ted llLorimenl200l| ). with or- 
bital periods as short as 3 hr JEd wards fc Bailea 2000'). The 
third class are white dwarf-black hole binaries. For black hole 
companions in the stellar mass range, there are no known 
examples-they would be extremely hard to find. There is, 
however, growing theoretical evidence for, and observational 
evidence consistent with, the existence of intermediate mass 
black holes (IMBHs) . X-ray observations of nearby galaxies 
hint at the possibility of a large population of black holes in 
the ~ 10^-10"^ Mr ma ss range JColbert fc Mushotzkvlll999l: 
IColbert fc Ptak^l2002^ ■ The interpretation of these observa- 
tions is the source of some debate, but accreting IMBHs are 
certainly a plausible explanation. These IMBHs may have 
been formed in the early universe from 10^-10^ M0 stars 
JBromm et alJll999r . Also, recent dynamical simulations of 
dense stellar clusters have demonstrated that runaway col- 
lisions can lead rapidly to the formation of black holes with 
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mass es in the ~ 10^-10^ M© range JPortegies Zwart et alJ 
|2004) . More generally, the population of IMBHs could span 
the entire interval between stellar and massive black holes 
(~ lO^-lO*^ M0). T hey could either tidally capture white 
dwar fs from the field JFabian et alJlQTStlPress fc Teukolskvl 
Il977h . or participate in three-body processes in stellar clus- 
ters and end up in bound orbits with white dwarf compan- 
ions teillJIlQT'i) . 

The tidal capture mechanism is simple in principle. 
When the radius of periastron is no more than a few times 
the Roche limit, the normal modes of the white dwarf can 
be excited non-resonantly. On egress, the energy of the 
oscillation is taken out of the orbit, which can become 
bound. The energy transfer need be no more than ~(30- 
100 km s-^f = 10""" erg g"^ = 10-100 eV per nucleon for 
capture to ensue. Subsequent periastron passages can lead 
to additional mode excitation and the possibility of heating 
the star if the oscillation energy is thermalised. For a black 
hole companion, this mechanism requires that the Roche 
limit lie not too far inside the event horizon, which sets an 
upper limit of ~ 10^ M© on the black hole mass. (When the 
black hole mass exceeds ~ 3 x 10^ Mq, it turns out that 
gravitational bremsstrahlung is more important than tidal 
capture and contributes a larger capture cross-section.) 

Whatever their detailed nature and origin, WDCO bi- 
naries evolve dynamically under the action of gravitational 
radiation. The orbital periods and eccentricities will de- 
crease until the former reaches the Roche period, ~ 10—100 s 
depending upon mass, when the star will be torn apart 
by tidal forces. During inspiral, a WDCO system will pass 
through a series of resonances between harmonics of the or- 
bital frequency and the white dwarf normal mode eigen- 
frequencies. Typically, the system will spend many orbits 
near each resonance, and consecutive resonances for a given 
mode will be separated by a much larger number of orbits 
associated with the gravitational inspiral time-scale. Passage 
through a sequence of such resonances will result in transfer 
of energy from the orbit to the oscillations and may drive 
the amplitudes of the oscillations non-linear, with possibly 
observable consequences. 

If the energy transfer from the orbit to the stellar os- 
cillation is as large as ~ 10 keV per nucleon and the mode 
energy is thermalised, then it is possible to detonate the 
star, releasing more nuclear energy (~ 1 MeV per nucleon) 
than the gravitational binding energy (~ 0.1-0.3 MeV per 
nucleon). Typically, for a black hole companion, the orbital 
binding energy will be greater than the nuclear energy re- 
leased, so most of the ejecta will be trapped. However, for 
a white dwarf or neutron star companion, the ejecta can 
become unbound. 

We note that WDCO binaries also c onstitute one of 
the prime potential sources for LISA (e.g. iFinn fc Thornd 
[2000). The considerations presented below may be relevant 
to gravitational wave observations of these systems. 

In order to understand how much energy is transferred 
from the orbit to the white dwarf during inspiral, it is neces- 
sary to consider the evolution of the system after it has lost 
much of its initial eccentricity. We present a general formal- 
ism for handling tidal resonances in WDCO binaries in the 
linear normal mode approximation. We restrict ourselves to 
non-rotating stars because white dwarfs are observed gen- 
erally to be slowly rotating, and they are not expected to 



maintain corotation during inspiral (c.f. iBildsten fc Cutler! 
Il992h . However, this needs further investigation. In addition, 
we mostly confine our attention to the (. = m = 2 /-mode 
because it is expected to be the dominantly excited one. 
Nonetheless, our formalism may be applied to other modes 
as well. Of particular interest may be g-modes, as these have 
lower frequencies than the /-modes, and can therefore be ex- 
cited at fundamental resonance in circular orbits before tidal 
disruption. 

We begin with a description of the general formalism 
in !J5| This includes a brief description of the resonant exci- 
tation of a simple harmonic oscillator (§2.1), a summary of 
the structure and normal modes of cold white dwarfs (§2.2), 
an overview of gravitational radiation reaction (§2.3), and 
our formalism for tidal excitation (§2.4). In §3, we discuss 
some physical aspects of the dynamical problem, including 
mode damping and time-scales. We also discuss the physical 
regimes where various effects are important. In §4, we apply 
the general formalism to calculate the energy transfer for a 
resonance in an eccentric orbit, ignoring back reaction. In 
§5, we discuss the regime where this result is expected to be 
valid, and describe the qualitative effects of back reaction. 
This is followed by a description of the long term evolution of 
WDCO binaries in the absence of back reaction. We present 
our conclusions in §6. 



2 GENERAL FORMALISM 

2.1 Simple harmonic oscillator 

Consider an undamped simple harmonic oscillator with nat- 
ural frequency wo and displacement x{t) subject to an ex- 
ternal force per unit mass F{t). The equation of motion. 



X + LJqX = F{t) , 

can be easily solved to get 

±{t) = K [C(t)] , 

x(f) = ^9 lat)] , 

where 

C(i)=Coe'"°'-fCiW, 
and 



Ci(t)EEe""°* / dt'e-"^"''F{t') 



(1) 

(2) 
(3) 

(4) 

(5) 



It follows from the expressions for x and x in terms of (^ that 
the total energy per unit mass of the oscillator as a function 
of time is given by 



E{t) = l\mf 



(6) 



Let the external force per unit mass now be of the form 

F{t)^Fo{t)cos[(t>{t)] , 

with the amplitude Fo and frequency (p being slowly vary- 
ing functions of time, and > 0. Resonance occurs when the 
relative phase of the driver and the oscillator becomes sta- 
tionary. This gives us the condition </>(t) — ujq- We assume 
that there is only one passage through resonance, and since 
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there will be little average energy transfer away from reso- 
nance, we need only consider the form of the forcing function 
near resonance. Let tn be the time when the resonance con- 
dition is satisfied, and expand the driver in a Taylor series 
around this point: 



F(tfl + r)~Fo(ti?)cos 



<j){tR) + ^or + (pitR) — 



(7) 



(Since the amplitude varies slowly with time, to lowest or- 
der, we can take the amplitude as constant through the res- 
onance.) With the definitions 



Fr = Fo(ifl) , (j)R = (^(ffl) , a : 



this becomes 



<P{tR) 



F{tR + t) ~ Fr cos 4'r + tjJoT + a 



2 2 



(8) 



The parameter a has the physical interpretation of being a 
measure of the fractional change in frequency over a char- 
acteristic period of oscillation. The requirement that the 
frequency of the driver is varying slowly therefore implies 
a <C 1. In other words, the driver can be considered har- 
monic with a well-defined frequency over several periods of 
the oscillator. We can also view a as a measure of the phase 
"drift" -i.e. a measure of how fast the driver accumulates 
additional phase. With this interpretation, it is easy to see 
that the time spent near resonance is given by (qljo)"^ , 
approximately. 

The problem of solving for the motion of the oscilla- 
tor is now mathematically similar to the theory of Fresnel 
diffraction, and the solution can be expressed in terms of 
Fresnel integrals. We relegate the detailed calculation to Ap- 
pendix^and simply summarise the result: 

Ci = -^e'^^o^.J Ji(T)exp(i(^H) 



-|-72(r)exp 



where I\ and I2 are given by expressions involving Fresnel 
integrals. The time-averaged energy per unit mass changes 
asymptotically by 



AE = 



4aujg 



1 + 2 



Eo 



TvFyiaiol 



cos I/; 



(9) 



where Eo — \(o\ /2 is the initial energy, and tp is an initial 
phase. Qualitatively, the velocity is in quadrature with the 
force well away from resonance, but the relative phase of the 
two becomes approximately stationary near resonance for a 
time interval ~ (qoiq)"^ , and there is a velocity change 
~ FR{aujo)~^''^. The presence of the ip dependent term re- 
flects the fact that the oscillator may gain or lose energy, 
depending upon its initial energy and the relative phasing 
with the driver near resonance. If we perform an ensemble 
average over initial phases (assuming a uniform distribu- 
tion), we find that the average energy transfer is given by 



(AE) 






(10) 



From Q it is then obvious that, for Eq <^ e, the initial 



phase is unimportant and the actual energy transfer will be 
very close to the average. The possibility of negative energy 
transfer only exists when 



Eo > 



4C2 



where C = cos tp. Note that, since the energy of the oscillator 
cannot be negative, it must be true that AE ^ —Eo- It can 
be shown that @ complies with this constraint. 

Simple, linear damping is conventionally treated by 
adding a term 27i to the left side of (0. When 7 ^ 
(aLLio)^ J the development of the oscillation will be uninflu- 
enced by damping, although the energy of the oscillation will 
be converted steadily into heat. However, when the damping 
is efi'ective on the time-scale of energy transfer, the ampli- 
tude of the oscillation will be reduced. Nonetheless, it can 
be shown that the energy that appears ultimately as heat is 
still given by B, independent of 7, as long as 7 ^ wo (see, 
for example, iLandau fc Lifshitall96g) . 



2.2 White dwarf oscillations 

In this paper, we confine our attention to homogeneous, 
non-rotating white dwarfs where the pressure is contributed 
solely by cold, degenerate electrons. Thermal corrections. 
Coulomb effects, as well as compositional discontinuities are 
ignored. The relevant equations of stella r structure are de- 
scribed in lKippenhahn fc Weigertilll99Cl) . We consider three 
cases with masses 0.6, 1.0, 1.4 Mq for fie = 2. Some relevant 
properties are given in Table 

The linear theory of normal modes for a cold 
white dwarf has been give n by several authors (e.g. 
iKippenhahn fc Weigertlll990l) . We briefly review it to es- 
tablish notation and to resolve some ambiguities. It is con- 
venient to use Eulerian perturbations and to choose four 
independent variables which we write as ri{r)Yf°^{6,if)e^'^'^, 
where 



r 



pgr 



gr 



(11) 



^P is the Newtonian potential, g = 'i/' is the gravitational 
acceleration, primes denote differentiation with respect to r, 
and Y/^ are the real spherical harmonics defined for m ^ 
by 



V° 
J- Im 



Ym , m = 

(Yi^ + Yl^) /V2 , m > 

, m = 

(y«™ - y'D /iV2 , m > 



Our choice to use real rather than complex spherical har- 
monics allows us to avoid complications that would arise 
from the expansion of a real vector field in terms of com- 
plex basis vectors. Note that the Cowling approximation, 
which ignores the potential variation, is commonly used to 
describe stellar pulsations, but it is inadequate for our calcu- 
lations. We denote the spatial displacement field for a mode 
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Mass 



Radius 



U)f2 



M 



/2 



(Mq) (10« cm) (w.) (10--^ M,) 



mj2{R,)/vij2(R*) 6/2 



0.6 


8.83 


1.53 


2.05 


1.0 


5.71 


1.65 


1.28 


1.4 


1.98 


1.97 


0.25 



-0.169 
-0.124 
-0.0412 



9.6 X 10^ 
8.0 X 10'5 
1.5 X 10^ 



Table 1. Homogeneous, cold white dwarf models with fie = 2, and properties of their quadrupolar /-modes. 



by ^j(r),^ the components of which are given by 



6(r) = 



'^'^'^1^ + Ji'^'^^^ ) ^^™(^' f) 



(12) 



With the above definitions, the continuity, Euler, en- 
tropy, and Poisson equation s yield four first-order ordinar y 
differential equations (e.g. iKippenhahn fc Weigerd flQQCh . 
There are two boundary conditions at the centre of the star, 



^1 = 



gV2 






and two more at the surface, 



Vi = V2 ~ V3 



V3 






With a choice of normalisation, the system of equations for 
Tf can be solved numerically to get the normal mode frequen- 
cies and eigenfunctions. The usual choice is to normalise the 
eigenfunctions by mass so that 



/ <f'xpo$.ri'j =Sj 



(see, for example, IPress fc Teukolskvl Il977f) . where po is 
the unperturbed mass density. However, we make the non- 
standard choice rii{R,) = 1, where R, is the radius of the 
white dwarf. This ensures that the coefficient of a mode in 
an eigenfunction expansion is a measure of the relative sur- 
face displacement and, consequently, a measure of the mode 
non-Unearity. With our normalisation convention, we define 
the mode "mass" by 

M, = ^J'd'xpom^ . (13) 

It is convenient for us to work in terms of dimensionless 
quantities whenever possible, and towards that end we shall 
make use of the definitions 



Rl 



P, 



GAU 
R, 



£* = 



GMj 



where M« is the mass of the white dwarf. The physical signif- 
icance of these definitions is that, to within factors of order 
unity, uj, is a typical mode frequency scale, /3, is the surface 
escape speed, and E^ is the gravitational binding energy of 
the white dwarf. 

^ A word on notation: we use the the index j as shorthand for 
the complete set of labels required to specify a mode uniquely. In 
places where we need to distinguish between even (oc cosimip)) 
and odd (oc sva(rrnp)) modes, we use the additional labels e and 
o for even and odd, respectively. In such places, j should be un- 
derstood to mean the set of all other labels. Occasionally, we use 
the index a which can take on the values e or o. 
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Figure 1. Radial eigenfunctions of quadrupolar /-modes for the 
0.6 M0 model from TableQ] 



The most important modes for our purpose are the 
quadrupolar /-modes. For a non-rotating star, the five /- 
modes with 1 = 2 are degenerate in frequency. The eigen- 
frequencies for our three white dwarf models are given in 
TableQ The radial eigenfunctions for the 0.6 Mq model are 
displayed in Figure^ The eigenfunctions for the other white 
dwarf models are qualitatively similar. 



2.3 Gravitational radiation 

We adopt a Newtonian approach to gravitational radiation 
reaction in the two-body problem, neglecting all finite-size 
effects. Namely, we treat the problem as essentially Keple- 
rian with prescribed corrections to the orbital equations. For 
non-relativistic orbits (v < 0.2c), the secular corrections due 
to gravitational radiation are provided to a fair approxima- 
tion by the orbit-averaged expressions 



dSo 



dt 
dLoTh 



32 ^ 

5 *'^*(l+g)2/3 



^Pt 



10/3 



^i(e) 



dt 



32^ q' 

5 (l-Hg)2/3 



pn- 



7/3 



^2(e) 



(14) 



(15) 



l|Peterall964|) . where i?orb and Lorb are the orbital energy 
and angular momentum, q is the ratio of the companion 
mass to the white dwarf mass, n is the Keplerian orbital 
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frequency, e is the orbital eccentricity, and 

We can re-express the orbital evolution in terms of changes 
in the orbital frequency and eccentricity: 

de 



5 *(l + g)i/3 



-LOt, 



V3 



—,P. [ — ] Tz{e) 



(16) 
(17) 



dt 15 (l + g)i/3 

where 

-r- / N e / 121 2 

■^3(e)=(^3^5)H7i(^l + ^^ 

If gravitational radiation is the only mechanism for orbital 
evolution, then it follows from these equations that 



e = ~g{e)- , 
n 

where 



(18) 



e(e) = 



19J^3(e) 



18JP"i(e) • 
Finally, we can integrate the above equation to get 



n(e) 



2\3/2 



n(0. 54101) 



gl8/19 



1 ^ 121 2 

1 H e 

304 



-1305/2299 



(19) 



As the orbit shrinks, it circularises, eventually according to 
e oc n^^ , approximately. 

For a more accurate treatment of gravitational radia- 
tion (especially for high eccentricities), and for the inclusion 
of other general relativistic effects, corrections to the or- 
bital acceleration can be added directly to the equations of 
motion. Detailed derivations and discussi ons of these correc- 
tions can be found in the literature (e.g. liver fc Will 19951) . 
and we shall not reproduce them here. 

It should be noted that it is not necessary to worry 
about relativistic apsidal precession as it will only affect ne- 
glected higher-order terms. 



2.4 Equations of motion 

Neglecting general-relativistic effects, the equations of mo- 
tion for a white dwarf-compact object binary can be ar- 
rived at via several approaches. A straightforward, Newto- 
nian "push-pull" approach to the excitation of tides and the 
resulting perturbations to the orbit is perhaps the most intu- 
itive. However, it is not transparently self-consistent in terms 
of conserving energy and angular momentum. A variational 
approach has the advantage of maintaining self-consistency 
explicitly, and a lso has the potential for exposing underly- 
ing s ymmetries IIGingold fc MonaghanlligSCI : iRathore et alJ 

liooi). 

We expand the physical displacement, ^, of fluid ele- 
ments within the white dwarf in terms of normal modes, 



C(r,t) = ^A,(t)$,(r) 



(20) 



E,^^M,Rl\A% 



\Aj^ + A,j„ + Uj [Aj^ + ^jo)J 



Note that, as mentioned in H2.2I we are not using the con- 
ventional normalisation for the normal modes. With our nor- 
malisation, the $,j have dimensions of length and the mode 
displacements Aj are dimensionless. The energy and angu- 
lar momentum in each even-odd mode pair can be written 
in terms of their displacements Aja (t) as 

1 

- — 1 

2 
and 

Lj = mMjRl (AjeAjo - AjoAj^^ 

JRathore et alJ l2003r , respectively. For an isolated mode 
(i.e. one not being excited), the following simple relation 
also holds: 

LOj 

The conserved energy and angular momentum for the entire 
system are given by 

E=-^[R +Ru)-—^- 

and 

L — fiR ii + 2_, Lj 1 



(21) 



(22) 



where R is the orbital separation, u is the angular coordinate 
in the plane of the orbit, Mq is the companion mass, fi is 
the reduced mass, and 



/,(R) = GMo /dVpol.-V 



1 



Ir-R| 



(23) 



jRathore et aUEooJ^. 

We can use 12111 as the Hamiltonian for the system, and 
with the definition Xj = R,Aj, it takes the simple form 

rr_PR , pI GMqM, 

2n 2^_R2 R 



+ E(^ + ^Af,-|., 



(A 

\2M, 



2 2 



^^ofj 



(24) 



where pR — /iR, pu 
menta conjugate to R, u, and Xj , respectively. We therefore 
see that the Hamiltonian is comprised of three pieces: the 
Keplerian terms for the orbit, a sum of independent har- 
monic oscillators for the normal modes, and a sum of terms 
of the form Xjfj which couple the modes and the orbit. The 
overlap integral fj therefore plays a dual role as a driving 
function for the excitation of tides, and in the disturbing 
function for the orbit. This is not surprising since the system 
is conservative. Hence, any energy and angular momentum 
transferred to tides must be extracted necessarily from the 
orbit. 

From Hamilton's equations for 1241 . we obtain the fol- 
lowing equations of motion: 



PR 



pi 



fiR^ 

x^ 
' du 



GMqM, 



dfi 
' dR ' 



(25) 
(26) 
(27) 



6 Yasser Rathore, Roger D. Blandford and Avery E. Broderick 



The terms involving the derivatives of fj give the perturba- 
tion of the orbit due to the excitation of tides, and we there- 
fore refer to them as the back reaction terms. The overlap 
integral fj can also be written as 



■'■'^ ' ^ '■> * \ R J I sm(mu) 

where 



(28) 



JRathore et all2003fl . and the bracket notation denotes that 
either cos(mu) or sin(mu) will be present. It will be useful 
for us to write fj in yet another way. From the usual Ke- 
plerian relation between the orbital frequency n and the 
semi-major axis a, it follows that 



/. 



qr]jM,i^lR^ 



, 2U+1)/Z ... 



(l + (7)(«+i)/3 \<jj,) \R) 

We now make use of the Fourier expansion 



cos(mM) 
sin(mu) 



-- j exp(imt;) = ^ Xf '' 

fc — — oo 



(e) exp{ikl) 



where v is the true anomaly, I is the mean anomaly (not to be 
confused with i), and the Fourier coefficients X^'™ (called 
Hansen coefficients; see Appendix 10 are real functions of 
the eccentricity. Noting that u = v + zu, where vj is the 
longitude of periapse, we have 



a \*'+i J cos(?nu) 
rJ \ sm{mu) 



E^. 



-(e+i-),' 



cos{kl + mzu) 
sin{kl + mzu) 

The overlap integral /, is therefore given by 
where 



(29) 



fc=0 



(30) 



■f]k - - 



2{^+l)/3 



q-qjM,ulRt, 
(l + q)(«+i)/3 



X^f. cos{kl + mzu) + X-f. coa{kl — mzu) 
XX sin{kl + mzu) — XZ. sin(fc/ — m,zu) 



(31) 



and we have used the shorthand X 



jk 



Xr-(t + l),m r 

X^j. , for econ- 



omy of notation. It should be understood in the expression 
for fjk that, for fc = 0, only the Xt, terms are present. 
For fc > 0, the Xrj_ terms can be combined using trigono- 
metric identities. However, it is simpler to note that, since 
X^'"^ oc e'*^"'"' , to lowest order in eccentricity, the X~f. terms 
will be suppressed by 2rn powers of eccentricity relative to 
the Xi terms. Therefore, for low to moderate eccentricities 
(< 0.6) and m > 0, the X~^, terms can be neglected to a 
good approximation. For the case m — 0, the X j. terms 
are identical. Hence, in all that follows, for m — one only 
needs to make the change Xt, — > 2Xt, . 

From the preceding discussion, we know that the driv- 
ing function fj for the excitation of a particular mode is 
an infinite sum of fjt terms. The phases that appear in the 
expression 13111 for fjk are all of the form kl ± m,zu. Thus, 



there exists the possibility of resonance whenever the rela- 
tive phase of the mode and one of these terms is stationary: 
Wj = kl i: mzu, where Wj is the phase of the mode. As 
mentioned previously, the kl — mzu terms will be suppressed 
by 2m powers of eccentricity relative to the kl + mzu terms. 
Thus, the dominant resonances will occur for Wj = kl + mzu. 
It might be thought that the above condition is equivalent 
to u!j = kn, but, in general, this is not the case. As the 
evolution of the orbit is dependent upon the tides via the 
back reaction terms in the equations of motion, there are 
complicated, non-linear dependencies implicit in each of the 
variables in the resonance condition. However, since we ex- 
pect the orbital corrections to be relatively small, it should 



be true that, at resonance, w. 



kn. 



3 PHYSICAL CONSIDERATIONS 

3.1 The tidal limit 

Clearly, our formalism for treating the evolution of a WDCO 
binary as a dynamical interaction between the orbit and the 
tides is only valid if the white dwarf is not tidally disrupted. 
In other words, we require that the white dwarf does not fill 
its Roche lobe. This requirement constrains the harmonics 
of the orbital frequency that a given mode can interact reso- 
nantly with. To quantify the constraint, we use the following 
approximation to the radius of the Roche lobe: 



rn 



0.49g 



-2/3 



R 0.6<j-2/3 -^ ln(l -h q-i/3) 
iEggletonlll983i) . It then follows that we require 



k> 



[0.6 + ?2/^ln(l + g-i/3)] 



3/2 



(l_e)3/2 



(l + g)i/2 



(32) 



where we have made use of the facts that the orbital separa- 
tion at periapse is a(l — e), and that tOj ~ kn at resonance. 
Note that we have implicitly assumed that the companion 
is more compact than the white dwarf, and hence is not 
disrupted. This is certainly true when the companion is a 
neutron star or a black hole. However, for the white dwarf- 
white dwarf case, the actual constraint is provided by the 
star that is disrupted first, which may be the companion. 

It should also be mentioned that the above approxima- 
tion for the radius of the Roche lobe is for circular, syn- 
chronous orbits. A more general treatment of the Roche 
problem may modify the tidal disruption regime. This is 
a possibility for future investigation. 



3.2 Importance of the £ 



2 /-mode 



The lowest £ modes that can be excited tidally are £ — 2. 
Modes with higher values of £ will have smaller overlap in- 
tegrals, since fj oc iZ"'^"*"^'. We may therefore infer that the 
primary modes that are excited outside the Roche limit are 
the £ — 2 modes. It is also the case that, with our choice of 
coordinates, the m = 1 modes will not be excited. This is 
easily seen by remembering that rjj oc P™(0), and 



PriO)^l ^ ' i£-m)\] 



£ + m even 
£ + m odd 
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(see, for example. lArfken fc Weberlll995r . Therefore, the 
only £ = 2 modes that are excited have m = 0, 2. Further- 

the m — modes will be sup- 



more, since X? 



Jk-m\ 



pressed by two powers of eccentricity relative to the m = 2 
modes. Hence, we deduce that the dominant modes for low 
to moderate eccentricities will have £ = m = 2. Also, since 
the p-mode frequencies increase monotonically with the ra- 
dial order, we can access (before tidal disruption) the lowest 
harmonic resonances for the modes with lowest radial order- 
the /-modes. 

Putting together the above considerations, we conclude 
that the mode excited with the largest amplitude in a cold 
white dwarf will be the £ = m = 2 /-mode. Note that in 
a warm star, ij-modes can also be excited. These will have 
lower frequencies than the /-modes, and their frequencies 
will decrease monotonically with the radial order. However, 
the structure of g-modes is sensitive to assumptions about 
the stellar model. If the modes are confined to surface layers, 
then the overlap integrals will be essentially zero, and the 
modes will not be excited tidally. 



3.3 Mode damping 

The formalism that we have presented in ^does not include 
any mode damping. In a realistic scenario, white dwarf oscil- 
lations will damp out over sufficiently long periods of time. 
While we shall mention some possible mechanisms through 
which this might occur, we make no attempt to provide an 
exhaustive analysis as there is an extensive literature that 
exists for this problem. 

Some possible mechanisms that have been considered 
for the damping of nonradial white dwarf oscillations include 
gravitational radiation, neutrino losses due to pycnonu- 
clear reactions, and radiative heat leakage (Osaki & Hansen 
|l973). The relative importance of each mechanism depends 
on the typ e of mode under consid eration, but it was demon- 
strated bv lOsaki fc HansenI il973r that the dominant damp- 
ing mechanism for quadrupolar /- and p-modes, in the linear 
regime, is gravitational radiation. However, their calculation 
contains a numerical error. We present a corrected deriva- 
tion in Appendix ini 

Another possible mechanism for the damping of modes 
with large amplitudes is by non-linear coupling to other 
modes. Thi s has been explored extensively in various con- 
texts (e.g. iDziembowski Il982l : iKumar fc GoodmanI Il996t 
IWu fc Goldreichll200 J) . and it has been shown that non- 
linear mode interactions can be important amplitude lim- 
iting effects. In this paper, we ignore this complication be- 
cause it is, in fact, one of our goals to study whether such 
non-linear amplitudes can be excited by passage through a 
sequence of tidal resonances in a WDCO binary. 

In stars with compositional discontinuities or solid in- 
teriors, turbulence may be excited at boundaries, which can 
lead to additional dissipation. 



orbit, is given by 

^ 5 {l + qf/-' 



UJt 



-8/3 



2560;* q 

iPeteral964l) . For an eccentric orbit, the inspiral time has to 
be calculated numerically. In general, for an eccentric orbit 
with a given period, this time is shorter by up to a factor 
of 1000 for eccentricities up to 0.9, relative to the circular 
orbit time. For eccentricities < 0.5, however, the circular 
orbit inspiral time is a fair approximation. 

The second relevant time-scale is the mode damping 
time. In general, the damping times for quadrupolar /- 
modes depend upon the white dwarf mass. Assuming grav- 
itational radiation as the mechanism, the damping time (as 
derived in Appendix [UJ is given by 

T, = ^/3-5^3-/(7?*) 



Mj , 4 



For our 0.6 Mq and 1.0 Mq models, this gives ~ 3000 
and ~ 100 years, respectively. Note that these are neces- 
sarily underestimates since our cold white dwarf models are 
highly centrally condensed. In contrast, the damping times 
for "moderately realistic" 0.4 Mq and 1.0 Mq models used 
bylOsaki fc Hansen (1973) are about 2.8 x lO'^ and 500 years, 
respectively. The damping times are therefore quite sensitive 
to the stellar model. 

Finally, the third time-scale of interest is the white 
dwarf cooling time. A rough estimate for this is provided 
by 

_ 4.7 X 10'^ years / M^/Mc ^ '^^"^ 



JKippenhahn fc Weigertr l990'l. where A is the atomic mass, 
and I/« is the white dwarf luminosity. For typical parameters, 
this gives a cooling time of ~ 10^ years, which is much longer 
than any other relevant time-scale. We can therefore ignore 
the thermal evolution of the white dwarf. 

In order for mode damping via gravitational radiation 
to be physically unimportant during the long-term evolution 
of a WDCO system, it is necessary that Tj > Tan- In other 
words, we require that the damping between resonances is 
negligible during the gravitational inspiral. This gives us the 
following constraint on the harmonics that we can consider 
for a particular mode: 

1 3/8 



fc< 



153671 



(l+<?)l/3 



13, 



^(^•) ( w. 



-4/3 



(33) 



where we have used the expressions for Tj and Tqr given 
above, and have also made use of ujj ~ kn at resonance. For 
our 0.6 A/q white dwarf model and mass ratios greater than 
a few, this constraint evaluates to 

.<iir^^^'^ 

Note that, for moderate to high eccentricities, this is overly 
restrictive, and the actual limit obtained from an evaluation 
of the inspiral time for eccentric orbits is higher. 



3.4 Time-scales 

For the long term evolution of a WDCO binary, there are 
several time-scales of interest to us. The first of these is the 
gravitational radiation inspiral time, which, for a circular 



4 RESONANT ENERGY TRANSFER 

Let us now consider a mode being excited resonantly on a 
white dwarf in an eccentric orbit around a compact compan- 
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ion. We shall neglect the back reaction terms in the equa- 
tions of motion, and hence the orbit can be taken to be Ke- 
plerian with corrections due to gravitational radiation (the 
validity of the no back reaction approximation will be dis- 
cussed in >I5.11 . We assume that we start exciting the mode 
resonantly at t = 0, and limit our analysis to the regime 
iit/n <^ 1, where n is given by 11611 . This is not particularly 
restrictive since the gravitational radiation time-scale is typ- 
ically much longer than the resonance time-scale. Finally, we 
shall also assume low to moderate eccentricities (~ — 0.5), 
and hence neglect the X~^ terms in 13111 . 

With the above assumptions, we can expand the or- 
bital elements and phases in 13111 in Taylor series around 
resonance (retaining only the zeroth-order term in the am- 
plitude) to obtain 

2(^+l)/3 



fjk - - 



qrijA'hu),R* 

(l+g)(« + l)/3 \^^, 



cos I 
sin ( 



i>jk + ojjt + kht^) 
bjk + i^jt + knt ) 



(34) 



where (pj^ is an initial phase. We now note that 1341 is ex- 
actly of the form of 0, with the identifications 



Fr^- 



qrijUj^Rf 



(l+g){f+l)/3 \M,J \lO, 



M, 



2(f+l)/3 



X 



jk > 



LOO 



2 



2kh 



(the division by Mj in Fr is necessary since it is fj/Mj that 
appears on the right hand side of ILJSl '). We can therefore 
immediately write down the resonant energy transfer: 



E, 



M. 



768 (1 + g)(2'?+i)/3 \M, 

(«-7)/3 (x 



jk 



kTl 



where we have averaged over initial phases. Using the fact 



-e,H. 



{AE,k) _ q 

E, (l + g)(2«+i)/3"'"' 

where the parameter 



fc e 



e, 



768 






(«-7)/3 



(35) 



(36) 



depends only upon the white dwarf model and the mode, 
and 

^_4(f-l)/3 



-jfe 



(e): 



:Fi 



i^tk) 



(37) 



contains all the dependence upon the eccentricity and the 
harmonic. The values of the parameter Oj for our 0.6 Mq, 
1.0 Mq, and 1.4 Mq white dwarf models are given in Ta- 
ble Q We see that the energy transfer decreases monotoni- 
cally (relative to the star's binding energy) with the mass. 
As Hjfc(e) ex e'^'*^"'"^, to lowest order in eccentricity, the 
energy transfer is typically a very sensitive function of the 
eccentricity. Also, for a circular orbit, it is clear that only 
the fundamental resonance, k = m, exists (as would be ex- 
pected on physical grounds). We remind the reader that, for 
m 7^ 0, the energy transfer given by 1351 is for a particular 
choice of even or odd mode. It should therefore be multiplied 
by a factor of two to obtain the total energy transfer to the 
even-odd mode pair. 



5 DISCUSSION 

5.1 Regime of validity 

We now consider in what regime, if any, the no back reaction 
approximation is valid. Qualitatively, we expect back reac- 
tion to change the orbital frequency as a mode is excited 
resonantly, which will tend to push the system away from 
resonance. Clearly, this will modulate the energy transfer 
at some level. However, if the change in orbital frequency 
is small compared to the resonance width, then we expect 
that the modulation of energy transfer will not be signifi- 
cant. On the other hand, if the change in orbital frequency 
is comparable to or larger than the resonance width, then 
back reaction will play a significant role. Another way of 
saying this is that the modulation of the energy transfer by 
back reaction is a second-order effect. Therefore, as long as 
the energy transfer is small enough, we are justified in ignor- 
ing back reaction. We can quantify this criterion by defining 
a resonance parameter 



Xjk-- 



An 



jk 



^K" 



(38) 



where Arijk is what the change in kn would be if the energy 
given by 1351 were to be taken out of the orbit, and Arf,^'^ is 
the resonance width. In general, we expect that for Xjk ^ 1 
back reaction will not play a significant role in modulating 
the energy transfer, where as for Xjk ^ 1 back reaction will 
be important. Using the estimate Ari^?^ ~ (2fcn)^", we find 



Xjk 



r- r- ,9-15/2 2 (8f-23)/6 

5 57r /3, %a] 

3 2048 gi/2(i_^g)(«-i)/6 

,_(«-7)/3 

i^tk) 



M, 

m7 



(39) 



Tt 



.3/2 



Figure |5| shows the numerical integration across a particular 
resonance for various values of Xife, both with and without 
back reaction. The first qualitative feature that stands out 
is that the energy transfer with back reaction tends to be 
smaller than that without back reaction. This is not surpris- 
ing since the system with back reaction is expected to spend 
less time near resonance. Quantitatively, we see that for this 
particular resonance with Xjfc ^ 0.1 we obtain nearly iden- 
tical numerical results with and without back reaction, with 
Xjk ~ 0.1 the results differ by a factor of order unity (about 
2), and with Xjk ~ 1 the energy transfers differ by an order 
of magnitude. 

The delineation of the back reaction and no back re- 
action regimes in the eccentricity-harmonic plane obtained 
with the above criterion for a quadrupolar /-mode on a 
0.6 Mq white dwarf and various companion masses is shown 
in Figure |3 It is seen that the region of parameter space 
where back reaction may be neglected, according to the Xjk 
criterion, grows with the companion mass. There is, how- 
ever, a reason to think that back reaction might actually 
play an important role in some regions of the parameter 
space where the Xjk criterion indicates otherwise. 

Consider the following thought experiment. Imagine 
that we are approaching a resonance with an initial phase 
that would lead to a net negative energy transfer in the no 
back reaction approximation. As we start removing energy 
from the mode and depositing it into the orbit, the orbital 
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Figure 2. The energy in the £ = m = 2 /-mode on a 0.6 Mq white dwarf is shown for a passage through the fc = 15 resonance 
with different values of the parameter Xjk obtained by varying the eccentricity, and with q = 10, 000. In each plot, the dashed curve 
is the system without back reaction, and the solid curve is the system with back reaction. The curves have been smoothed to remove 
high-frequency components. 



frequency will necessarily decrease (i.e. the semi-major axis 
will increase) and the system will get pushed off resonance. 
It will then have another chance to approach the same reso- 
nance. Then, if the phase is such that energy is transferred 
to the mode, then the system will once again get pushed 
off resonance, but this time in the opposite direction (since 
the orbital frequency will increase). Gravitational radiation 
will then evolve the system away from this resonance and 
towards the next one. This scenario hints at the possibility 
that back reaction may force the resonant energy transfer 
to be always positive. However, this is not necessarily the 
case. For instance, we have assumed that there is sufficient 
initial energy in the mode to be able to change the orbital 
frequency significantly. Also, we have neglected the fact that 
gravitational radiation will be removing energy from the or- 
bit as we are transferring energy to the orbit from the mode. 
If the rate of dissipation by gravitational radiation is high 
enough, then back reaction may not matter. The system will 
evolve through resonance regardless, on a time-scale deter- 
mined by the rate of dissipation. Hence, we can still get a 
net negative energy transfer to the mode. 

In summary, back reaction may be important in deter- 
mining both the magnitude and the direction of resonant 
energy transfer. The Xjk criterion provides, in some sense, 
only a measure of the correction to the magnitude. In the 



regime where Xjk ^ 1, the implication is unambiguous: back 
reaction will be essential in determining the energy transfer. 
However, when Xjk < 1, things are somewhat uncertain for 
reasons stated above. A solution of the problem including 
back reaction is required to determine conclusively whether 
back reaction is important in that regime. 



5.2 Long term evolution 

In we calculated the energy transfer for an individual 
resonance in the absence of back reaction. In general, as the 
binary shrinks under gravitational radiation, the system will 
pass through a sequence of resonances for each mode. How- 
ever, this is only a possibility for an eccentric orbit because, 
as demonstrated previously, only the fundamental resonance 
exists for a circular orbit. We note that, in the no back reac- 
tion approximation, the energy transfer at a resonance can 
be negative as well as positive, depending on the relative 
phase of the mode and the driver, and the initial ampli- 
tude. Also, there will be negligible average energy transfer 
between resonances, as long as we are well outside the tidal 
limit. If we assume (as seems reasonable) that the system 
has no long term phase memory, then the relative phasing 
at each resonance will be essentially random, with a uniform 
distribution. It then follows that, on average, the mode will 
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Figure 3. The regions in eccentricity-harmonic space where back reaction is and is not important (labelled as "BR" and "No BR", 
respectively) are delineated according to the Xjfc criterion for a <? = m = 2 /-mode of a 0.6 Mq white dwarf, and various mass ratios. In 
each plot, the solid curve traces out the contour Xjk = li and the long dashed lines to its left and right trace Xjfe = 0.1 and Xjfe = lOi 
respectively. The short dashed lines trace three gravitational radiation inspiral trajectories through the plane. For reference, the tidal 
limit and the region where mode damping via gravitational radiation during inspiral is important arc also shown. Note that, for the 
mode damping curves, the full, eccentricity-dependent inspiral times have been used, rather than the circular orbit times. 
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tend to gain energy over time, and that the average total 
energy transfer after a sequence of resonances will be sim- 
ply the sum of the individual average energy transfers given 
by (El. 

Let £k denote the average energy transfer given by II35II 
for a particular mode at the fc-th resonance, and let Ek be 
the energy (in units of B*) in the mode before the fc-th 
resonance. It then follows from @ and our assumptions of 
random phases and negligible energy transfer between res- 
onances that, for a sequence of resonances in the no back 
reaction approximation, the evolution of the mode energy 
will be given by the discrete random walk (with a drift) 



Ek-i — Ek + Ek 1 + 2 




(40) 



where Ck is a random variable drawn from the distribution 
1 



p(x) = 



2^6 [-1,1] 



7r\/l — x'^ 

For a derivation of elementary statistical properties of this 
random walk, see Appendix IHl 

Figure |1] shows the results from calculations of passage 
through a sequence of resonances performed using the above 
random walk model for several sets of initial conditions. We 
have chosen to plot the mode amplitude 



B, 



1 

it* 



2Ej 



(41) 



rather than the energy, because we want to draw attention 
to the fact that, for moderate initial eccentricities, the am- 
plitude of a ^ = 7n = 2 /-mode can be driven to values in 
the range ~0.1-1. (An amplitude of unity for a £ = m = 2 
mode corresponds to a maximum physical displacement of 
the stellar surface of about 55% relative to the unperturbed 
radius.) We therefore expect that the linear normal mode 
analysis might not be valid in those cases, and that non- 
linear effects may in fact determine the actual outcome. 

It should be noted that, even if back reaction plays a role 
in determining the direction of energy transfer, our result 
that non-linear amplitudes for a ^ = m = 2 /-mode can 
be attained by passage through a sequence of resonances 
is unlikely to be affected. This is due to the fact that the 
result depends chiefly upon the allowed magnitude of energy 
transfer, and as we restricted our calculations to the regime 
where Xjk <Si 1, back reaction is not expected to change 
things. 



6 CONCLUSIONS 

In this paper, we have discussed a specific dynamical prob- 
lem: what is the energy transfer to the normal modes of a 
white dwarf due to tidal resonances in a compact object bi- 
nary inspiralling under gravitational radiation? For simplic- 
ity, we have considered only non-rotating, completely de- 
generate stellar models. We have provided a semi-analytical 
answer for the case when the perturbation of the orbit by 
the modes may be neglected (the no back reaction approxi- 
mation). It has been shown that, ignoring back reaction, and 
assuming that the relative phase of the mode and the orbit 
is randomised between successive resonances, the energy in 
a mode executes a discrete random walk, with a drift, during 



passage through a sequence of resonances. In addition, we 
have found that the amplitude of a ^ = m = 2 /-mode can 
be driven into the non-linear regime as a result of passage 
through such a sequence. We have also demonstrated that 
back reaction can modulate the amount of energy trans- 
ferred significantly, depending upon the parameters of the 
system. Our attempt to delineate the regime in parameter 
space where the no back reaction approximation is valid has 
had mixed success: we can predict when back reaction will 
be definitely important, but a clear delineation of the regime 
where it will be unimportant has not been found. Despite 
this ambiguity, our result that non-linear amplitudes can be 
attained for a. £ = m = 2 /-mode is expected to be robust. 
A treatment of the problem including back reaction will be 
the subject of a future publication. 



ACKNOWLEDGMENTS 

RDB would like to thank J. P. Ostriker, K. S. Thorne, 
P. Goldreich, Y. Wu, and M. L. Aizenmann for useful discus- 
sions and encouragement. The authors acknowledge support 
under NASA grant 5-12032. 



REFERENCES 

Abramowitz M., Stegun I. A., eds, 1972, Handbook of 

Mathematical Functions. Dover, New York 
Arfken G. B., Weber H. J., 1995, Mathematical Methods 

for Physicists. Academic Press, San Diego 
Bildsten L., Cutler C, 1992, ApJ, 400, 175 
Bromm V., Coppi P. S., Larson R. B., 1999, ApJ, 527, L5 
Colbert E. J. M., Mushotzky R. F., 1999, ApJ, 519, 89 
Colbert E. J. M., Ptak A. F., 2002, ApJS, 143, 25 
Dziembowski W., 1982, Acta Astronomica, 32, 147 
Edwards R., Bailes M., 2000, ApJ, 547, L37 
Eggleton P. P., 1983, ApJ, 268, 368 
Fabian A. C, Pringle J. E., Rees M. J., 1975, MNRAS, 

172, 15P 
Finn L. S., Thorne K. S., 2000, Phys. Rev. D, 62, 124021 
Gingold R. A., Monaghan J. J., 1980, MNRAS, 191, 897 
Hills J. G., 1977, AJ, 82, 626 

Iyer B. R., WiU C. M., 1995, Phys. Rev. D, 52, 6882 
Kippenhahn R., Weigert A., 1990, Stellar Structure and 

Evolution. Springer- Verlag, Berlin 
Kumar P., Goodman J., 1996, ApJ, 466, 946 
Landau L. D., Lifshitz E. M., 1969, Mechanics. Pergamon 

Press, Oxford 
Lorimer D. R., 2001, Living Reviews in Relativity, 4, 5 
Maxted P. F. L., Marsh T. R., Moran C. K. J., 2002, MN- 
RAS, 332, 745 
Misner C. W., Thorne K. S., Wheeler J. A., 1973, Gravi- 
tation. W.H. Freeman and Co., San Francisco 
Murray C. D., Dermott S. F., 1999, Solar system dynamics. 

Cambridge Univ. Press, Cambridge 
Osaki Y., Hansen C. J., 1973, ApJ, 185, 277 
Peters P. C, 1964, Physical Review, 136, 1224 
Portegies Zwart S. F., Baumgardt H., Hut P., Makino J., 

McMillan S. L. W., 2004, Nature, 428, 724 
Press W. H., Teukolsky S. A., 1977, ApJ, 213, 183 



12 Yasser Rathore, Roger D. Blandford and Avery E. Broderick 



2 


- ' ' 










1 1 1 


M 




■"x 


q=10 


3 


- 




'\ 


- 








^^ 


- 


4 






^\, 


- 


5 


- 






- 


6 


- 










Ill 


1 - 





















1 








s 


q=100 




-2 


_ 


\ 




e„=0.1 


- 








\ 








^^ 






^ 








\ 




-4 




\ 






" 


m" 


- 




"■^ 


o 






^^ 




- 












- 






_ 


n'"- 


- 




-6 




^ 










^^ 












^^^ 


- 






1 1 1 1 1 1 1 1 1 1 1 


, , 1 , , 





8 10 12 14 16 

k 



6 B 10 12 14 IB 

k 



q-1000 
e„=0.2 









^ 














1 


1.5 


- 




'^ 


q=1000 
e„=0.3 










^^ 


^ 


-2 


- 












- 










'\ 










2.5 


- 








- 


-3 


- , 














1 " 



10 



14 16 18 20 

k 



q=100,000 
e„=0.2 




q-100.000 
e„ = 0.4 



Figure 4. The amplitude of a £ = tti = 2 /-mode of a 0.6 Mq white dwarf during passage through a sequence of resonances in the no 
back reaction approximation is shown for several sets of initial conditions. All of these lie in regions of the eccentricity-harmonic plane 
where back reaction is not important according to the Xjk criterion. The calculations were done using our semi-analytical formalism. 
In each case, the solid line shows a particular realization of the random walk given by 1401 . and the dashed line follows the ensemble 
average. The random walks were terminated when Xjk ~ 0.01. Note that the scales on the axes are different for each plot. 
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(see, for example. lAbramowitz fc StegurJll972f) . Let us now 
consider the energy after long times. If we start out from 
zero frequency, then luotq — —a, and 

1 . 1 



<-l 



> 1 



APPENDIX A: DETAILS OF CALCULATION 
FOR FORCED HARMONIC OSCILLATOR 



From JSJ and Q, it follows that 



Ci = ^e dr { exp 



+ exp 



i\ (t>R + a 






i (jifl + 2ujqt' + a 



2 II 



Completing the square in the argument of the second expo- 
nential, we get 



Ci = —e"^''''{ Ii{T)exp{i(j)R) 



+ /2(r)cxp 



^ \ (pR 



where 



hi'^) = / d'^' exp ( -uloIt'^ 



hir) = / dr' exp 



U r- '^2 

2 \ Ja 



Making the change of variables 

u{t') = J —loot' , 

I I. _ JOL , 2 

V TT VTra 



we have 



du exp 



12 = — \ — "■'" exp 

^0 V a Jvo 

where uo — u(to), Ut = u(r), et cetera. Evaluating the 
integrals gives 



The I2 integral will not contribute significantly since its up- 
per and lower limits are effectively the same. For 7i, we find 

1 /tt ,, . .^ 1 /27r_,;^/4 



/i = — ,/-(l+J) = — \/— e 
Hence, after some algebra, 



AE = E- Eq 



AolOq 



1 + 2 



Eo 



nFl/Aaul 



cos ip 



where -Bo = lCo| /2 and -0 is an initial phase. 



APPENDIX B: HANSEN COEFFICIENTS 

The Hansen coefficients X^'™ for the two-body problem are 
defined by 

RY 



— j exp(imw) = y^ X^''"(e) exp(ifc/) , 



where R is the orbital separation, a is the semi-major axis, v 
is the true anomaly, e is the orbital eccentricity, and I is the 
mean anomaly. The Hansen coefficients are real functions of 
the eccentricity, and it can be shown that, to lowest order 
in eccentricity, 



JMurrav fc Dermottll999l . and references therein). The coef- 
ficients can be calculated to any desired order in eccentricity 
as a series in terms of Newcomb operators; 



XI'- [e) = e 



I fc — m I 






+c 



where A = max(0, fc — m), C, — max(0,m — k), and the 
Ne wcomb operators X°', are d efined via recursion relations 
fsee iMurrav fc Dermotti 1999). Alternatively, for quantita- 
tive work, the Hansen coefficients can be evaluated for a 
given eccentricity by calculating the integral 

^r„fRv 

2-K 



Xl'^^ie) = ^ I dl {-] cos{mv - kl) 



numerically. 



h 



1 /tt 



(jjo V '^ 
1 17 



{CK) - C{uo) + i [S{Ur) - S{uo)] } 



h = —\l-\CiVr) - Civo) - i [S{V^) - Sivo)] I 

where the Fresnel integrals are defined by 



C{z) 



S{z) 



Note that 



dt cos 



dt sin 



TTi^ 



■nr 



lim C(z) = lim S{z) = ±- 

2 — vzboo z — »itoo Z 



APPENDIX C: DAMPING OF QUADRUPOLAR 
MODES UNDER GRAVITATIONAL 
RADIATION 

The average power radiated in gravitational waves due to 
a time-dependent mass quadrupole moment is given, in the 
weak-field limit of general relativity, by 

dt 45c5 ^^'^^'^^ 

JMisner et al.lll973l) . where Qij is the mass quadrupole mo- 
ment as defined conventionally in classical physics: 



Q^. 



jd\ [,..., -r^&.,)pi^,t) 
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It can be shown that the power radiated by a quadrupo- 
lar mode is independent of m (this is a consequence of 
the Wigner-Eckart theorem). Hence, we can restrict our- 
selves to the m = case for simplicity. It then follows 
that Qii — Q22 ~ —(533/2, and that the off-diagonal terms 
vanish. Therefore, noting that the time dependence is sinu- 
soidal, we have 



dK 



GW 



dt 



G 
' 60c5 ' 



.6>S2 



(CI) 



where Q33 should now be understood to mean the time- 
independent amplitude of the mass quadrupole moment. 
Writing the mass density as 

p{x) = po{r) + BjSpj{r)Y2o{r) , 

where Bj is the amplitude of the mode and Spj is the nor- 
malised density perturbation associated with the mode, we 
find 



Q33=4J-B, 



-R. 



drr 5p{r) 



The above integral can be simplified by using the linearised 
Poisson equation, integrating by parts twice, and using the 
surface boundary condition ri4 = —{£ + l)i]3. The result is 



^33 = \j-M,Rlri3]{R*)Bj 



Substituting the above expression into ICIL we get 



-4,{R^J^Bi 



dEX^ GMiRt 

Finally, noting that the total energy for an isolated mode is 
given by 

Ej = Ui.RW.B^, , 



Since E^ will only depend upon Ca, for a > k, Zk is linear 
in Ck. As all the Ck are independent random variables (by 
assumption), it follows that 



{Zk)^0, {Zl)=2ek{Ek) 



(D2) 



Given an initial mode energy Ek before the fe-th res- 
onance, we wish to determine the average mode energy 
(Ek-p), and its variance ap, after passage through p res- 
onances. From (IDIL we have 



Ek-p = Ek + ^ Sa + Ap , 
a — fc — p-t-l 

where we have defined 

k 



(D3) 



ip _ 7 ^ Za 

a^k — p-\-l 



Using <D2ll . it follows immediately that 

fe 
(Ek-p) ^Ek+ J2 ^P- 

/3^k-p+l 



(D4) 



Note that, since e/3 > 0, (Ek-p) increases monotonically as 
we pass through a sequence of resonances. Hence, we say 
that the random walk JDll has a drift. To calculate the 
variance, we need to find (Ek-p)- Writing 




we see that 

{El_p) = {Ek-pf + {AD , 



(D5) 



as all the terms linear in Ap will vanish when averaged. To 
calculate (Ap), we write 



we arrive at the e-folding time for the mode energy under 
damping by gravitational radiation: 



Gtt 



/^.-V/(i^.)(^)^. 



J 



Y, zi + 2 Y. E ^"^' 



13 



d^fc— p+1 



a = fe-p+l I3 = a+1 



and note that, since /3 > a in the above double sum, each 
term of the double sum will be linear in Ca , and will, hence, 
vanish upon averaging. Therefore, we have 



APPENDIX D: STATISTICAL PROPERTIES OF 
RESONANT ENERGY TRANSFER IN THE NO 
BACK REACTION APPROXIMATION 

In this appendix, we derive some statistical properties of the 
random walk given by 14UII . which we rewrite as 

Ek-i = Ek+ek+Zk , (Dl) 

where 

Zk = lyskEkCk ■ 

Recall that Ek is known in advance, and that Ck is a random 
variable drawn from the distribution 



-pix) = 



ttVT 



, xe[-l,l] 



Note that 



(a) = o, kcl) = - 



a— fc— p+l 

Substituting into 1D5L and then using IID4II . we find 

k / '^ \ 

al = 2 J2 £c\Ek+ J2 ^A ■ 



(D6) 



a^fc— p+l 



/3 = a+l 



It should be noted that lIDlll is not Gaussian, nor will 
it become Gaussian after many resonances. That the pro- 
cess is not Gaussian is clear from the fact that the random 
walk is bounded from below. Furthermore, the central limit 
theorem is not applicable because, typically, the probability 
distribution of Ap is dominated by the most recent few har- 
monics, and hence the effective number of variables never 
becomes large. 

This paper has been typeset from a Ti^X/ I^TijiK file prepared 
by the author. 



